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Hard-Loop Dynamics of Non-Abelian Plasma Instabilities 
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Non-Abelian plasma instabilities may be responsible for the fast apparent quark-gluon thermaliza- 
tion in relativistic heavy-ion collisions if their exponential growth is not hindered by nonlinearities. 
We study numerically the real-time evolution of instabilities in an anisotropic non-Abelian plasma 
with an SU(2) gauge group in the hard-loop approximation. We find exponential growth of non- 
Abelian plasma instabilities both in the linear and in the strongly nonlinear regime, with only a 
brief phase of subexponential behavior in between. 



In this Letter we present first results on the real- 
time evolution of non-Abelian plasma instabilities due 
to momentum-space anisotropies in the underlying quark 
and gluon distribution functions P, S S S '"i 7] in the 
nonlinear hard-loop approximation. Such anisotropics 
are generated during the natural expansion of the mat- 
ter created during a heavy-ion collision and the result- 
ing instabilities may be responsible for the fast appar- 
ent thermalization 8], which seems to be faster than 
can be accounted for by perturbative scattering processes 
0, 0, 01 . This type of plasma instabihty is the 
analogue of the electromagnetic Weibel instability which 
causes soft gauge (magnetic) fields to become nonper- 
turbatively large. Eventually this leads to large-angle 
scattering of hard particles 13], thereby rapidly acceler- 
ating the isotropization and subsequent thermalization of 
an Abelian plasma with a temperature anisotropy. How- 
ever, in the non-Abelian case it is conceivable that the in- 
trinsic nonlinearities could cause the instabilities to stop 
growing before they have large effects on hard particles 
and therefore reduce their efficacy in isotropizing a non- 
Abelian plasma. 

The regime where the backreaction of collective fields 
on the hard particles is still weak but where the self- 
interaction of the former may already be strongly nonlin- 
ear is governed by a "hard-loop" effective action which 
has been derived in Ref. 6] for arbitrary momentum- 
space anisotropies ^20j- We discretize this effective ac- 
tion in a local auxiliary-field formulation, keeping its full 
dynamical nonlinearity and nonlocality. This is then ap- 
plied to initial conditions that allow for an effectively 
l-|-l-dimensional lattice simulation, extending a previous 
numerical study 7] that used a static and linear approx- 
imation to the hard-loop effective action. 

Discretized Hard-Loop Dynamics. — At weak gauge 
coupling g, there is a separation of scales in hard mo- 
menta IpI — of (ultrarelativistic) plasma constituents, 
and soft momenta ~ g|p| pertaining to collective dynam- 
ics. The effective field theory for the soft modes that is 
generated by integrating out the hard plasma modes at 
one-loop order and in the approximation that the ampli- 
tudes of the soft gauge fields obey <^ |p|/ff is that 



of gauge-covariant coUisionless Boltzmann-Vlasov equa- 
tions |14|. In equilibrium, the corresponding (nonlocal) 
effective action is the so-called hard-thermal-loop effec- 
tive action Is'l which has a simple generalization to plas- 
mas with anisotropic momentum distributions Its 
contribution to the effective field equations of soft modes 
is an induced current of the form 
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W^{x;^), (1) 



where / is a weighted sum of the quark and gluon distri- 
bution functions The quantities Wp{x]v) satisfy 



[v ■ D{A)]Wp{x-^) ^ Fp^{A)v^ 



(2) 



with = p'^/IpI = (l,v), = d^- ig[A^,-], met- 
ric convention (H ), and this has to be solved self- 

consistently with D^{A)F^^'^ = j'^ . At the expense of 
having introduced a continuous set of auxiliary fields 
Wp{x;v) the effective field equations are local, but non- 
linear in the case of a non-Abelian gauge theory. 

We shall be interested in particular in the case where 
there is just one direction of momentum-space anisotropy 
(e.g., the collision axis in heavy- ion experiments), so we 
assume cylindrical symmetry and parametrize the first 
derivatives of the hard particle distribution function by 
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Because p^Wjs = 0, we have 
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v^'[fiW°ix;v)+ hW'ix-v)]. (4) 



In the isotropic case one has /2 = 0, and only W'^ 
appears, whose equation of motion |(5J) is driven by 
(chromo)electric fields; in the anisotropic case, however, 
enters, whose equation of motion involves the z com- 
ponent of the Lorentz force. 
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The induced current QJ is covariantly conserved as 
one can verify by partial integration with respect to 
0, but no partial integration is needed for parity in- 
variant distribution functions /. This is only a mild 
restriction on the choice of /'s, but a great advantage 
for the following discretization. For /(p) = fi^p), 
the two terms in the current Q are covariantly con- 
served individually, because Eq. JSjl implies that D ■ j (x 
J d^p {fiF^^v^ + f2F^''v^), which vanishes by symmetry. 

Our aim will be to study the hard-loop dynamics in 
lattice discretization where we approximate the contin- 
uous set of fields W/3(a;;v) by a finite number of fields. 
Because Eq. Q does not mix different v's, a closed set 
of equations is obtained when the integral in Eq. is 
discretized with respect to directions v 

^ E^^'t^v^W + 6vMC(x)], (5) 

V 

where the TV unit vectors v define a partition of the unit 
sphere in patches of equal area, and where , 6v are then 
fixed coefficients for a given distribution function /. Co- 
variant conservation of j is ensured by symmetric choices 
of the set of v's which satisfy J2v = 0, = 0, 

and J2v bv^± — 0, with v^ = v — v^e^. 

Given such a discretization, which is a discretization 
of the phase space of the hard particles with respect to 
the directions of their momenta, only Af auxiliary fields 
Wv = a^W^ + bvW^ participate in the dynamical evolu- 
tion. The full hard-loop dynamics is then approximated 
by the following set of matrix-valued equations, 

[v ■ D{A)]W^ - {a^F'"' + b^F'^')v^ (6) 

D,{A)F>^- = j^Y.''-W^^ (7) 

V 

which can be systematically improved by increasing A/". 

A different possibility for discretizing hard-loop dy- 
namics that has been employed previously in isotropic 
plasmas would be a decomposition of the auxiliary fields 
W^{x\ v) into spherical harmonics p^H^. Our proposal 
is slightly simpler, but also more flexible in that it allows 
one to e.g. improve approximations in highly anisotropic 
cases with cylindrical symmetry by selectively increasing 
the resolution in z direction more than in Lp direction. 

The dynamical system described by Eqs. 10 and {Tj) 
has constant total energy of the form 

J d^x tiiE"^ +B^)\t+ dt' J d^x 2 ti it' ■'Ef- (8) 
The part containing the induced current involves 



which shows that in the isotropic case (/2 = 0) there 
is a local, pos itive definite energy contribution from 
the plasma [l9j- However, in the general anisotropic 
case, positivity is lost, corresponding to the possibility 
of plasma instabilities, where energy may be extracted 
from hard particles and deposited into the soft collective 
fields without bound as long as the hard-loop approxi- 
mation <g; \p\/g remains valid. 

1+1-dimensional lattice simulation. — When only z- 
dependent initial conditions are imposed, the entire dy- 
namics proceeds l-|-l-dimensionally and we can take all 
collective fields as 1-Fl-dimensional, though the underly- 
ing hard degrees of freedom are of course still 3-1-1 dimen- 
sional, with a discrete set of directions v for their (con- 
served) momenta. This dimensionally reduced situation 
already allows us to study the evolution of non-Abelian 
Weibel instabilities, which in the linear (Abelian) regime 
are formed by (superpositions of) transverse standing 
waves with exponentially growing amplitudes. 

We have solved the equations © and Q in this di- 
mensionally reduced situation by lattice discretization, 
closely following Ref. who considered a toy model cor- 
responding to an induced current which is simply 

j^^^fi^A^, a^x,y. (10) 

As was shown in Ref. 7] , this reproduces the static limit 
of the anisotropic hard-loop effective action for fields that 
vary only in the anisotropy direction, but it neglects its 
general frequency dependence and dynamical nonlinear- 
ity. This is now provided through Eqs. © and Q. 

Like Ref. we work in temporal axial gauge A'^ = 
and take initial conditions corresponding to small ran- 
dom chromoelectric fields E = —dtA with polariza- 
tion transverse to the z-axis, and all other fields vanish- 
ing, which in our case includes the auxiliary fields Wv 
This initial condition satisfies the Gauss law constraint, 
D • E = A/""^ yVv5 whose continued fulfilment is mon- 
itored in the simulation, but not enforced. As a further 
nontrivial check of our numerics, we tracked conservation 
of the total energy ((SJ. 

In order to be able to compare with the analytical re- 
sults of Ref. , we consider an anisotropic distribution 
function / = TV(0/iso(p^ + CPz)- This determines the 
coefficients in Eq. (O according to 

av = TO^(1 +C^'^)"^ b^^^v^a^, (11) 

where = A^(^)m|,(0) with muiO) the Debye mass 
of the isotropic case ^ = 0, and N{^) a normalization 
factor which equals N{^) = ^/l + ^ if one requires that 
the number density of hard particles remains the same 
for different values 
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FIG. 1: The growth rate 7(fc) of unstable modes for the dis- 
cretized hard loops with M — 100 (dashed line) in comparison 
with the continuum result (solid line) for ^ = 10. 



For ^ > transverse modes with wave vector |k| < n, 



u? — - \ 1 + ^ — arctan \/l 
4 V 



(12) 



are unstable and grow with a fc-dependent growth rate 
7(fc), shown in Fig. ^ for ^ = 10, which has a maximum 
7* at nonvanishing |k| — fc* and a zero at k = 0. By 
contrast, the growth rate implied by the toy model H10() 
is given by 7 = -^//i^ — k-^, whose maximum is at k = 
and equals /i. 

In the linearized case one can also determine the dis- 
persion laws analytically in the case of finite M . Fig. 
compares the growth rates 7(fc) for full and discretized 
hard loops with J\f ^ Uz ><■ M'ip = 20 x 5 = 100, which 
shows that the latter give very accurate approximations 
for M > 100. As another example, the asymptotic mass 
of the stable propagating modes with |k| ^ is given 
for full hard loops by — m? /[2-/S^) arctan v^, which 
for ^ = 10, TV = 100 is reproduced with an error of only 
0.017%. 

The results of our lattice simulations of the nonlin- 
ear evolution for a gauge group SU(2) are finally shown 
in Figs. 121 and |21 where we have used the parameters 
^ = 10, Nz X A/'cp = 20 X 5 for discretizing the unit sphere 
with uniform spacing in z and Lp. The one-dimensional 
spatial lattice has 10,000 sites, periodic boundary condi- 
tions, and lattice spacing corresponding to anioa = 0.028 
so that the physical size is L « 280 m'^. Using a leapfrog 
algorithm with time steps e/a = 1/100, we track the evo- 
lution of a single field configuration |2lj with random 
initial seed chromoelectric field of root-mean-square am- 
plitude 0.012 m^l g. The Gauss law constraint turns out 
to remain preserved within machine accuracy; the total 
energy JSJ is conserved within less than 1%, with the 
percentage level being reached only at the largest times 
when the energy in the soft fields has grown by more than 
a factor of 10^ (full details will be given elsewhere). 

The upper panel of Fig. [3 shows the evolution of 
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FIG. 2: The average current jVms and the average relative size 
C of commutators [jx , jy] as a function of time (upper panel); 
the Abelianization correlation length ^aIJ] and the ordinary 
correlation length ^^[j], normalized to the scale of maximal 
growth A* (lower panel). 



which (after some initial wobble) grows exponentially 
with a growth rate that is most of the time only slightly 
below 7* , except for a transitory reduction at the begin- 
ning of the nonlinear regime, when jrms ^ "^iL/ff- Also 
shown is the dimcnsionless observable C, defined by |2^ 



C 



L tr (j2 + j2) 



1/2 



(14) 



(13) 



and giving a measure of local "non-Abelianness" . In the 
toy model of Ref. it was found that C suddenly be- 
gins to decay exponentially when fields get strong, with 
a decay rate similar to the growth rate of jrms- In the 
hard-loop case, we observe a similar phenomenon, but the 
decay rate of C is much smaller (and also less constant 
once the decay begins). 

Ref. also observed a concurrent global Abelianiza- 
tion, by comparing the correlation among parallel trans- 
ported spatially separated commutators to the correla- 
tion of parallel transported fields. A correlation length 
S,A defined through the former was found to rapidly grow 
to lattice size when C begins to decay, whereas no such 
growth occurred in the general field correlation length 
In the lower panel of Fig.|21we show the analogous quan- 
tities computed using currents instead of field values [22] 
and we found that correlated Abelianization takes place 
over extended domains, which remain bounded, however. 
In fact, turns out to be comparable with the scale of 
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FIG. 3: Average energy densities E in transverse/longitudinal 
chromomagnetic/electric fields and the total energy density 
contributed by liard particles, f (HL). 

maximal growth A* = 27r/A:* w Irri^ . By contrast, in 
the model of Ref. 7] A:* vanishes, which is presumably 
responsible for the different global behavior. 

Fig. 121 shows how the exponentially growing en- 
ergy transferred from hard to soft scales is distributed 
among chromomagnetic and chromoelectric collective 
fields. The dominant contribution is in transverse mag- 
netic fields, and it grows roughly with the maximum rate 
7* both in the linear as well as in the highly nonlinear 
regime, with a transitory slowdown in between. Trans- 
verse electric fields behave similarly, and are suppressed 
by a factor of the order of (7*/^*)^ |23. The appearance 
of longitudinal contributions, which are absent in the ini- 
tial conditions we have chosen, is a purely non-Abelian 
effect. While completely negligible at first, they have a 
growth rate which is double the one in the transverse 
sector, and they begin to catch up with the latter just 
when local Abelianization sets in. At that stage there 
appears to be some complicated rearrangement taking 
place, which delays the exponential evolution by a time 
Ai ~ 7~^, but subsequently the growth rate gets restored 
roughly to its previous value. Eventually there will be a 
point where the hard-loop approximation breaks down, 
namely when the energy transferred from hard to soft 
modes becomes comparable to that initially present in 
the former. 

From the above numerical results on the hard-loop dy- 
namics of non-Abelian plasma instabilities we conclude 
that the latter do not saturate until they begin to have 
large effects on hard particle trajectories, even though 
our simulations indicate complicated dynamics and only 
limited Abelianization of the unstable modes. Therefore 
it appears indeed possible that non-Abelian plasma in- 
stabilities are responsible for accelerated thermalization 
in a weakly coupled quark-gluon plasma. 
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